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Abstract 

In the present paper we compare the theory of mixtures based on Rational 
Thermomechanics with the one obtained by Hamilton principle. We prove that 
the two theories coincide in the adiabatic case when the action is constructed with 
the intrinsic Lagrangian. In the complete thermodynamical case we show that 
we have also coincidence in the case of low temperature when the second sound 
phenomena arises for superfluid Helium and crystals. 

1 Introduction 

The first mathematical model of homogeneous mixture of fluids in the context of Rational 
Thermodynamics was due to Truesdell [lj. The compatibility with the second prin- 
ciple of thermodynamics was well established by Muller in the framework of classical 
mechanics [2] and by Hutter and Muller in relativity [3] . 

In the framework of binary mixture of Euler fluids, Dreyer [H [5] was able to revisit 
the well known Landau model of superfluidity [HI [7j. The second sound phenomena in 
the case of liquid He II is now well explained from a macroscopic point of view. Recently 
Ruggeri [8] observed that a mixture of two Euler fluids can be regarded as a single heat 
conducting fluid. This result is advantageous to explain the second sound phenomena of 
crystals with the same model than for superfluid helium. 

A different approach was given by Gavrilyuk et al, [9], Gavrilyuk and Gouin [TQl 
[TT] . They consider a variational approach to describe two- velocity effects in homogeneous 
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mixtures: a Lagrangian of the system is chosen as a difference of the kinetic energy of the 
two constituents and a volumic potential which is Galilean invariant depending on the 
relative velocity of components. The equation of motions of the two components are not 
in balance form (in fact they are in balance form in Lagrangian variables associated with 
each component). Nevertheless, the momentum and the energy equations for the total 
mixture are in the clasical balance form. 

The present work compares the previous approaches and proves that the two theories 
coincide in the mechanical case when the Hamiltonian action is constructed with the 
intrinsic Lagrangian, i.e. does not depend on the relative velocity. Such is the case with 
the Lagrangian considered by Gouin in [12J. In the thermodynamical case we prove 
also the coincidence in the case of low temperature and we obtain a complete agreement 
between the two approaches and the superfluid model considered first by Landau. 



2 The Binary Mixtures of Euler Fluids 

The thermodynamics of a homogeneous mixture of n constituents is well codified as a 
branch of Extended Thermodynamics |13j . It is based on the metaphysical principles 
of Truesdell [1] which postulates the same balance laws of a single fluid for simple 
mixtures. 



2.1 The Balance System 

The equations of balance of mass, momentum and energy of the constituents read as 
follows 

+div(p a V a ) = T a , 



dt 
dt 



+ div (p a v a (giv a - t ) = m a , (a = 1,2, ...n), (1) 



d(\p a vl+ Pa^a) f/1 \ 1 

— Qj. + div | y-p a v a + p a e a j v a - t a v a + q a | = e a . 

These equations have the same form as the balance equations for a single body, except 
for the non-zero right hand sides which represent the production of masses, momenta and 
energies. These productions are due to interaction between the different constituents. Of 
course, since the total mass, momentum and energy of the total mixture is conserved, we 
must have 

n n n 

a=l a=l a=l 

where p a , v OJ e a , t tt , q a are the mass density, velocity, internal energy, stress and heat flux 
respectively of the a-component of the mixture. 
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If we sum the equations ([I]) over all constituents and introduce 

n n 

the density p = ^2p a , the velocity v = — v a , (2) 

0=1 a=l P 

the diffusion velocity u a = v a — v, (3) 

n 

the stress tensor t = (t a — p a u a ®u a ) , (4) 

a=l 

n 

the intrinsic energy density psj = p a £ a , (5) 

a=l 

1 n 

the internal energy density pe = pej H — Pa -1 ^ (6) 



" f 1 1 

and the heat flux q = ^ \ q a + p a (e a + -u 2 a )u a - t a u a V , 

a=l ^ 1 > 



(7) 



we obtain for the total mixture: 
The balance mass 



^ + div (pv) = 0, 



The balance equation of momentum 

dpv 



dt 

The balance of energy 



div (p v ® v - t) = 0, (9) 



d(lpv 2 + pe] r/i \ i 

V ^ + div|f-pv 2 + pej v-tv + q| = 0. (10) 

Note that equations ([HI [U [TOl) have the same form as those for a single fluid. Moreover in 
equation (1101) for the balance of energy we observe that the total kinetic energy is -^pv 2 
is not the sum of the kinetic energy of the components. In fact we have 

1 1 n 1 n 

a=l a=l 

By analogy with the intrinsic internal energy we call intrinsic kinetic energy the expression 

1 n 

Z a=l 

As we consider a single absolute temperature T, the aim of extended thermodynamics for 
fluid mixtures is the determination of the 4n + 1 fields : 

mass densities p a 

velocities v a (a = 1, 2, ... n). 

temperature T 
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To determinate these fields we need an appropriate number of equations. They are based 
on the equations for each constituent of balance of mass momentum (pQ)2 and con- 
servation of energy of the total mixture ( FlOl) . 



2.2 The Equations of Binary Mixture of Euler Fluids 

We consider a binary mixture of Euler fluids, i.e. fluids that are neither viscous nor 
heat-conducting : 

q a = 0, t a = -p a I, (a = 1,2). 

Instead of the mass and momentum balance laws for the second component, we use the 
equivalent equations of total conservation for mass and momentum. Therefore, associated 
with the 9 unknown fields (p 1; p 2 , Vi, v 2 , T), we have the 9 balance equations: 

^ + div (pv) = 

dpi i. / \ 

+ dlV (PiVi) = T\ 



Ot 
<9pv 

~dT 

dpjvi 
dt 



+ div(pv<g>v- t) = (11) 
div (PiVidgVi + pil) = ni! 



d(\pv 2 + pe) r/i \ -J 

— i — — '- + div | i-pv 2 + pe\ v - tv + qj = 



with 



q = E {a. ( £ » + ^) + p*} u «' 

2 

t = - E (pJ+P a u a ®u a ) , (12) 

a=l 

2 



a=l 



2.3 The Entropy Principle and Thermodynamical Restrictions 

The compatibility between the system (TT]) and the entropy principle expresses in the form 

^ + div {pSv + > 0, (13) 
or 
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which yields several restrictions on the constitutive equations [13] 



pS = Pl Sx + p 2 S 2 (14) 
Pi = Pi(pi,T);p 2 =P2(p2,T);e 1 =e 1 (p 1 ,T);e 2 = e 1 (p 2 ,T) (15) 



such that 



TdSi = de\ — ^\ dp x ] TdS 2 = de 2 — ^§ dp 2 (16) 
Pi P2 

q 1 

* = T ~ f ^ Pl/ilUl + ^2 U 2) • (17) 

where p a = e a + — — TS a is the chemical potential of constituent a. 

Pa 

2.4 The Mixture considered as a Single Heat conducting Fluid 

Ruggeri [8] proved that it is possible to write the velocities of the two constituents in 
terms of mass velocity and heat flux centers : 

a a 
Vi=vH q, v 2 = v q 

Pi P2 

where 

a V Pi 2 / V P2 2 / 

Introducing the concentration c = — , equations (IHp ? and fll 11) ^ can be written in terms 

P 

of p, c, v and q and the system (1111) becomes: 

2jL + div (pv) = 

°(P C ) , s 

— h div (pcv+aq) = r 

dpw ( a 2 \ 

-^- + div pv<g>v +pl+— - -q®q =0 (19) 

at \ pc(l — c) J 

d(pcv-\-o.c\] f c? I 

h div < pcv ® vH q ® q+a (v(g>q + q(g>v) + z/I> = — bq 

at [ pc J 

d(\pv 2 + pe) fn x / a 2 v x i 

1 drv<M-pv 2 + p£ + p v+ — + l q =0. 



at [V2 r r v ^(i-c 
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To eliminate the index 1, we write as in |8i], v — p\,r — T\ and mi = — £>q. In an 
extended thermodynamic model with 9 fields, the binary mixture can be considered as a 
single heat conducting fluid with a variable concentration. 

Equation of evolution (fl9]) /i is a natural extension of the Cattaneo equation for the 
heat flux. Thermal inertia term a together with term v have to be interpreted as new 
constitutive functions. The advantage of this procedure comes from the fact that the two 
functions are now understandable in the light of mixture theory: term v plays the role of 
one-component pressure while the thermal inertia term a given in (fl8l) is the inverse of 
the difference between the non-equilibrium enthalpies of the two constituents. 



2.5 The Superfluidity and Second Sound 

Dreyer [4J proved that the Landau theory of superfluidity is a particular case of simple 
mixtures with the thermodynamical peculiarities : 

S s = 0; p s - p n + ^(v s - v n ) 2 = 0, m s = r s v s , (20) 

where the indexes n and s correspond to normal and the superfluid components. 
By neglecting the quadratic term in the second equation, in the small diffusion case the 
two chemical potential p s and p n must be equal. Consequently, the relation p s = p n 
allows to obtain one field variable in terms of the others and it is possible to write 

Ps = P S (P,T) 

In this case equation ffTTj^ evaluates the mass production value r s and the superfluid 
helium framework becomes a theory with 8 fields (i.e. the system is formed by equations 
(HI])i,(III])3,(III])4,(in])5 or equivalents equations PI) i,® 3, ©4, (111]) 5 )• 

The condition (TSUjU is the most complex. In fact ([TTj) 4 with (fTTT) ? can be rewritten 
(see [S] for details) : 

dv s „ ( 1 9 \ 

-gf + V i^-v 2 s + fi s j + curl v s x v s = 0. 

This equation is in balance form only when the involutive constraint curl v s = holds. 
In this case the system f lT9"j) coincides with the Landau model [6J : 

^ + div(pv) = 0, 



dpv . 
-^j- + div (pv <g> v — tj 

<9v, _/l 



(21) 



+ V ( ~-v 2 s + p 



dt \2 



d(lpv 2 + pe) m r n 



dt 



+ div | {i^pv 2 + pej v - tv + q| = 0. 
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Taking into account fl24l) . ffTTl) and fl20|) i , the entropy law reduces to the Clausius form: 

fipS 'div f^v+lUa (22) 



dt V T 

where the heat flux (|T21) i is: 

q = pTS u n + i (p s u 2 s u s +p n u 2 n u n ) . (23) 

In the diffusion velocity we neglect the third order terms and we obtain the Landau 
entropy law for the heat flux [6]. The entropy flux becomes pS v n and the entropy is 
convected by the normal component 

^ + div (pSv n ) = 0. (24) 
at 

To focus on the thermal wave associated with the second sound we consider a rigid body 
at rest with constant density. For the superfluid component, the system of energy and 
momentum equations is : 

^ + div q =0, 

9Vi +vfk+^=o, 



dt V2 



with q = pTS v n . Such a system is in the form (fT9|) for a single fluid : 

dps 

— + chvq=0 
dt 

The system coincides with the one deduced by Ruggeri and coworkers for the model 
of second sound in crystals [14J. Such a model explains the change of form of the initial 
square thermal waves both in crystals [HI [H3, [16] and in the superfluid helium |17j . 



3 The Hamiltonian Procedure for Two-Fluid Mix- 
tures 

To obtain the equations of motion and energy, the procedure is the following: 
Let us suppose that the mixture of two miscible fluids is well described by the two- 
component velocities vi, v 2 , the densities p x , p 2 and the intrinsic internal energy [3 = psj. 
The intrinsic internal energy is a Galilean invariant and does not depend on the reference 
frame. We consider the general case where (3 depends on p l1 p 2 but also of the relative 
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velocity w = vi — V2 through the norm u — |vi — V2I [9]. The intrinsic kinetic energy is 
E c = 2 (pi^i + ^2) • 

Without dissipative effects, chemical reactions and with conservation of masses of the two 
components, an extended form of Hamilton principle of least action is used in the form 

51 = with I — L dxdt, 
JWo 

where the Lagrangian is L = E c — (3{p 1 , p 2 ,u), W = [t , ti] x D is a time-space cylinder 
and the variations must vanish on the boundary of W. The virtual motions of the mixture 
are defined in [9], [TO] . 

From the variations of Hamilton action, we obtain the equations of motions in the form 
^+curlk a xv B + V^-i«2 + k B v a ^ =0 (a =1,2) (25) 

where 

k . = v ._(_i ) .±WE. 

p a du u 

The momentum conservation law is obtained by summing on a = 1,2 equation (1251) 
multiplied by p a : 

( <9/3w£*Dw\ 
+ div p 1 v 1 (g) vi + p 2 v 2 ® v 2 - - = (26) 

\ ~ OLO UJ J 

d(3 

Additive terms come from the dependance of (3 in u and in the mechanical case p-,— h 

dpi 

9(3 n 1 
p 2 -^ p represents the total pressure p. 

dp 2 

The conservation of energy is obtained by summing on a = 1,2 equation ( 1261) multiplied 
by p a v a : 

^ + 2^ + ^ + "fa) + dlV l, PlVl ^ + klVl + P ^ 2 dp-2 + J = ( 7) 

In paragraph 2, we consider the case where /3 is independent of u and the entropy principle 
(TT5|) presented in [T3] yields /? = p 1 £i(/9 1 ) + p^ip^)- Then, equation (123]) writes 

+ curl v a x v a + V (- v 2 a + /i a j =0, (a = 1, 2). (28) 

Multiplying equation (1281) by p a straightforward calculations yield equation (IITL with 
m a = 0. Equations (1261 1271) yield equations flTTlk . (ITTjk and balance of mass equations 
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correspond to r a = (a = 1, 2). 

A purely mechanical case is the adiabatic one and we have verified the following results : 
In the adiabatic case with intrinsic Lagrangian L = E c — pei difference between the in- 
trinsic kinetic energy and the intrinsic internal energy pej = Pi^i(pi) + P2 e i{.P2)i the 
system deduced from Hamilton principle coincides with the system coming from Rational 
Thermomechanics. 



4 The Hamiltonian Procedure for Superfluid Helium 

In the case of a binary mixture some change must be done in the definition of virtual 
motions presented by Serrin in [TSj. Let us consider the motion of Helium II as two 
diffeomorphisms 

z = M(Z), z = M n (Z n ) 



where z 



corresponds to the Eulerian variables in time-space and Z 



,Z n 



correspond to the Lagrangian variables associated with the barycentric motion 



and the normal component motion of helium II. In coordinate form, 

9 ( A , X ) A /l/T {7 \ _ ( 9 (A n , X n ) 



M(Z 



M n ZJ 



4>(x,x)J , "" " \(p„ i\„. x„; 

We consider three one-parameter families of virtual motions which are sufficient to obtain 
the governing equations : 

( t = g (A,X) = g n (A n ,X n ) 
(?) x=d>(A,X, £ ) 

I x = <t>n (An,X n ) 

with $(A,X,0) = 0(A,X), 

{t = p(A,X) = g n (X n ,X n ) 
x = 0(A,X) 
x = (A n ,X n ,e) 

with $ n (A n , X n , 0) = (p n (A„, X„) , 

!t = G (X, X, e) = G n (A n , X n , e) 
x = 0(A,X) 
X = Cp n (A n , X n ) 

with G (A, X, 0) = G n (A n , X n , Q)=g (A, X) = g n (A n , X n ). 
The three families generate the virtual displacements 







( ° 1 



£=0 





An. 



( \ 



> St 



Ct 



e=0 



/dG\ 

de 

V o ) 



6 = 
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The virtual motion (J 7 ) generates an associated displacement 5Z n of the normal compo- 
nent. Indeed, the relations 

g(X,X) = g n (X ni X 1 ) 
e/> n (A,X n ) = $(A,X, £ ) 



imply 



c 



f dgn dg n \ 
d(f> n d(f) n 



V 



5Z r . 



J 



By using the definition of the deformation gradient F proposed in Appendix we get 



/ 



5Z n — C n £ 



with 











\ r n v n 5 r n 



(29) 



In the same way, virtual motion (jF n ) generates an associated displacement 5 n Z of the 
barycentric motion 



SnZ = CC n 



with 





F X V 





F-i 



Now, H (Z, e) notes a perturbation of h (Z), the variation of h is 

dH 



5h 



e=0 



We can also introduce Lagrangian variations corresponding to the families {J r n ) and {Tt 



de 



and 



dH t 



e=0 



de 



e=0 



The variations of the entropy S is a main step of our model: we make the physical 
assumption that the entropy S is defined on the Z n -space. This result corresponds to 
equation (l2"lj) proposed by Landau. Consequently, we deduce S n S = and 5 t S = 0. 
From relation fT2"9~|) we obtain 

cn dS OS 

6S= dz n 6Zn= d^t 

Following the Hamiltonian procedure presented in paragraph 3, we consider the La- 
grangian L as a function of p,v, p n ,v n , S (L = L(p, v, p n , v n , S)). Such is the case 

for the intrinsic Lagrangian L = —{p n v^ + p s v 2 s ) — (3(p,p n , S) where p s and v s are given 
by the relations : 



Ps = P - Pn 



and 



PV ~ Pn^r, 
P- Pn 



(30) 
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Consequently, 

<9v s 1 dv s 1 dv s p dv s p n 

-«- = — (v-v s ), — = — (v a -v n ), — = — I, — = I 

dp p s dp n p s <9v Ps dv n p s 

The variation of the Hamilton action corresponding to the first family is : 

51= j 5(L det B) dw 

JWq 

dz 

where B = —— is the Jacobian of M and Wo is the associated Lagrangian domain in 



the ^^.^j -space. Consequently, 



51= \ (5L + L Div C) detB dw . 

JWo 



Variations of L come from 



with, 



rr dL. dL . 8L £ 8L . dL :a 
5L = — dv + - — dv„ + — 5p + -w~5p n + wt;5S. 
dv dv n op dp n dS 

dL dL . . 

^=PV S , _ = p n (v B -V.) 
f)T 1 

R = Tp = -ifs + v s v - P' Pa ( Pn , p., S), (31) 

Rn= ^~ = + \ V s - VsVn ~ {3' Pn (p n , P s , S) + f3' ps (p n , Ps , S), 

pT = ~ds 

Moreover we have, 

x dv n dv n dp n dp n 

s ^ = d2r n CnC = ^t and 5p « = dz n CnC = ^t 

Since £ — (^^j ' we ( see Appendix for the variations 5p and 5v variations), 

n , T n . , dLd£ dL dv n dL 

5L + L Div C = -TT-jT + ~ 7T- € ~ P -7T div £ 

(TV at av„ ax op 

^ dL dp n ^ + L dw ^ + dL dS 
dp n dx. dS dx 

d£ , \ d\ n „ „ . . . „ dp„ ,. r . . . dL dS . 
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By using the expression 

pv s — = — (pv s (P v s) £ + div(p(v ® v s ) £) - div (pv ® v s ) £ 

we get 

5L + L div £ = 

^ (pv s £)-q£ (pv s ) £ + div(p(v ® v s ) £) - div (pv ® v s ) | 
+ P„v n ^^- div(pi ^ ) + V (pi?)£ + i? n ^ + div(L£) +|§V5| 

and from equations fl3TT) . 

5L + L Div C = 
(-J^(pv s ) - div (pv ® v s ) + V ( P R) -RVp-p v)j £ 

d 

+ 7^7 (pv s C) + div(p(v ® v s ) t) - div (pR£)+ div (L |) , 



at 

where * notes the transposition. Consequently, the first equation of momentum is 

^ + V [-v 2 s + (3' ps j = v s x curl v s (32) 
If we note p s = (3' p , when v ~ , equation f[3"2"j) yields 

at 



+ V 7-^ + p s =0 (33) 



which is the Landau equation for the superfluid component. In fact Landau pointed 
out that Helium II lose its superfluidity when the velocity is not small enough and the 
supplementary term curl v s x v « corresponds to this experimental evidence. 

Variations of the Hamilton action are closely the same for the second family. The 
variation of the entropy is 5 n S = and consequently an entropy term is now appearing 
in the equations of motion. The second equation of momentum is 

d 

T^(Pn( V " - V s)) + dlV ( Pn V n ® (V n - V 8 )) 

+ Pn(^) (v„ - v s ) - Pn VR n - p T VS = (34) 
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By summing equations fl32|) and (1311) . equation flM|) can be replaced by the balance of 
total momentum: 

^(pv s + p n (v n - v s )) + div (pv <g> v s + p s v n g) (v„ - v a ) - p|^ - p n ^- + Lj = 0. 

Straightforward calculations yield the equation of momentum 

<9pv 

-7^- + div (pv n ® v n + pv s (g) v s + p) = 0, (35) 
where p = p s p s + p n p n — (3 is the total pressure, with p n = j3' p . 

Finally, the third family is associated with the vector displacement <^ t = ( J . The 



variations of basic variables are calculated in Appendix : 



£ t v = -v^, 5 t p = pVrv, 5 t v n = -v n ^j^, 5tp n = p n VTv n , 8 t S = 0. 

The variation of the Hamilton action is 

5 t I = [ ( 5 t L + L^- ) det S rfw 
Vw \ at / 



with 



Hence, 



. T dL . aX . aX . aX . aX . 
<>X = ^-c) t v + - — 5tv n + -^-dtp + T—o t p n + -^-o t s 
av av„ op op n os 



, T rdr fdr „ \ . , f Or ^ \ 

b t L + L— = -pv s v I — + Vr v I - p n (v n - v s )v n I — + Vr v n 1 

+pRVr v + p nJ R„Vr w n + — (Lr) - — r 

at at 

= ~di ^ VsV ^ + di ^ VsV ^ T _ div (p( v s v ) v r ) + div(p(v s v) v)r 
'^(Pn( v n - v s )v„r) + ^(Pn( v n - v s )v n )r - div(p n v n (v n - v s )v n r) 
+ div p n ((v n - v s )v n )v n t + div (pRv t) - div (pi? v) r 
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+ div (p n R n v n T) - div (p n R n w n ) t + — (Xr) - — r. 

Consequently, 

d 

^ (pv s v + p n (v n - v s )v„ - X) + 
13 



div {(v s v - R) pv + [(v n - v s )v„ - R n ] p n v„} = 

If we notice that 

pv s v + p n {y n - v>„ — L — -p n v 2 n + ^p s v 2 s +(3 = pe 

and 

(v s v - R)pv + ((v n - v s )v n - f?„)p n v n = 

we obtain the equation of balance of the total energy in the form : 

dps 



dt 



+ div q = 0. (36) 



We notice that the specific entropy S does not appear explicitly anymore in equations (I32I) . 
(I35p . (1361) and we conclude : in i/ie case o/ superfluid Helium the Hamilton principle 
yields the LANDAU model. 



5 Appendix. Variation of Basic Tensorial Quantities 

Let (A, X) be any generalized Lagrangian coordinates and (£, x) the associated Eulerian 
coordinates 

x = 0(A,X). {67} 

The relation dx = v dt + FdX defines simultaneously the velocity vector and the defor- 
mation gradient of motion (j3"7|) : 

dd) 1 dd> d(f) dg 1 

v ^ 



d\ dg_ ' ax 8\dXdg_' 

OX d\ 

Let { ^ ^"L. , »' "^r \ be a virtual motion. The associated perturbation of the velocity 
lx = *(A,X,e) 

v is given by the formula : 

_ d<f> 1 

U " ~dX dG 

OX 

and consequently, 



s ^ U 

de 



d£ 1 <9r 1 

£=0 ~ dx~dl~ Y dxW 



dX OX 

For fixed values of Lagrangian coordinates the variation of v in Eulerian coordinates is : 

d£ dr d d d 

ov = — v— where — = — + v— . 

dt dt dt ot ox 
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Analogous calculation for F is : 



5F 



( 




Moreover, the Euler-Jacobi identity yields 



b det F = det F tr (f^f) . 



Hence, the mass conservation law is : p det F = p (X) and implies 



bp = —p ( div £ — Vr ■ v) 



Equation (jTTj) ? is the form of the mass balance for the normal component of Helium. If 
we assume 

p n det F n = p 0n (A„,X n ), 

which means that p n is defined on the Lagrangian space of the normal component, the 
variation of p n with respect to b n is always in the form : 
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